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Abstract 

The focus in this paper is on elliptic homogenization of a certain kind 
of possibly non-periodic problems. A non-periodic and two-dimensional 
example is studied, where we numerically illustrate the homogenized ma- 
trix. 

1 Introduction 

Background. When studying the microscale behavior (beyond the reach of 
numerical solution methods) of physical systems, one is naturally lead to the 
concept of homogenization, i.e., the theory of the convergence of sequences of 
partial differential equations. 

The homogenization of periodic structures using the two-scale convergence 
technique is well-established due to the pioneering work by Gabriel Nguet- 
seng [To] and the further development work by Gregoire Allaire [1]. Gener- 
alizations of the two-scale convergence technique have been developed indepen- 
dently by, e.g., Maria Luisa Mascarenhas and Anca-Maria Toader (scale 
convergence) , Gabriel Nguetseng [TTl |T2] (S-convergence) , and Anders Holm- 
bom, Jeanette Silfver, Nils Svanstedt and Niklas Wellander [31[B] ("generalized" 
two-scale convergence). 

A simple but possibly powerful method of analyzing non-periodic struc- 
tures is the A-scale convergence technique introduced by Anders Holnibom and 
Jeanette Silfver [5]. A-scale convergence is scale convergence in the special case 
of using the Lebesgue (i.e.. A) measure and test functions periodic in the second 
argument 5 . Homogenization techniques based on this approach are developed 
in the doctoral thesis jTS] of Jeanette Silfver. These results are the point of 
departure for the main contributions in this paper. 



Organization of the paper. In Section [2] we look at the convergence for 
sequences of functions. We start by stating the definition of the traditional 
notion of two-scale convergence as introduced by Gabriel Nguetseng [IlJj . Then 
we move on to generalizations of this convergence mode, namely "generalized" 
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two-scale convergence [H |S] and scale convergence [5] . We conclude the section 
by introducing A-scale convergence and the important notion of asymptotically 
uniformly distributed sequences [51 |T3] . 

Section [3] deals with the convergence for sequences of partial derivatives. We 
first look at how the two-scale convergence works for partial derivatives in the 
periodic case, and then we consider the more general case of A-scale convergence 
of sequences of partial derivatives. 

The concept of homogenization is introduced in Section 21 which begins 
with stating the definition of H-convergence [S], i.e., the generalization of Sergio 
Spagnolo's G-convergence of sequences of symmetric matrices [T31[TS]. We give a 
theorem on the homogenization of a sequence of periodic matrices. We then in- 
troduce the important "type-H^" property, which is employed in the end of the 
section when formulating a theorem on the homogenization of A-structures |13j . 

In Section [5] we specifically study a non-periodic and two-dimensional exam- 
ple of a A-structure with the property of not only having a periodic direction, 
but also having oscillations with a monotonically increasing frequency in one, 
non-periodic direction. We formulate and prove a series of proposition which 
are needed in order to prove the main result of this paper, Theorem|30l in which 
we present a homogenization result for the considered A-structure in the form 
of a homogenized matrix and a governing local problem. 

The concluding Section [HI illustrates the results achieved in Section [5] We 
numerically solve the local problem to obtain the non-constant homogenized 
matrix, and we heuristically explain why it is isotropic on a line along the 
periodic direction. 

Notations. The following more or less handy notations are employed: 

The A^-tuple (^i, . . . , ^n) in is denoted {^i)i=i,...,N , or whenever conve- 
nient. Similarly, (mij)^ jv, or simply the mere majuscule M when handy, 
denotes an TV x TV matrix. A bold dot, i.e. • , represents a non-fixed variable, 
e.g., (f>{- ,y) is the same as the function x i—s- 4'{x, y) where y is held fixed as a 
parameter. We will also allow expressions like, e.g., k • meaning x i~> kx. In this 
paper, O is always an open bounded non-empty subset of and, if nothing 
else is stated, Y is the unit cube (0, 1)^ in R^. The TV-tuple {■^X=i n 
partial derivative operators is denoted by the symbol V. Following, e.g., [3], 
the function space M^pcr(^) denotes the subspace of functions in H^^^{Y) with 
vanishing mean value. 

2 Convergence for sequences of functions 

The two-scale convergence method was introduced in 1989 by Gabriel Nguet- 
seng |10) , and a modern formulation is given by Definition [T] [7j . 

Definition 1. A sequence {u'*} in L'^{Q) is said to two-scale converge to the 
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limit uq e L^(ri X Y) if, for any v G Cpor(i^)), 

lim j u^{x)v{x,hx)dx ~ J J UQ{x,y)v{x,y) dy dx. 
n n Y 

An important property of the two-scale convergence is given by Proposi- 
tion [2 [2]. 

Proposition 2. If {u^} two-scale converges to uq, and v £ L^{^', Cpci{Y)) , 
then 



u'* uoi-,y)dy inL^{n) 

Y 



and 

vi-,h-)^ J v{-,y)dy inL\n). 

Y 

Analogous to a corresponding compactness result for weak convergence, we 
have Theorem [3] [7]. 

Theorem 3. Every hounded sequence in L^(ri) has a subsequence which two- 
scale converges. 

Introducing sequences of operators defined below we may extend Defi- 
nition [T] of two-scale convergence to a generalized version according to Defini- 
tionHHIS]. 

Definition 4. Assume that Y is an open bounded subset of M*^. Let X C 
L^(fl xY) be a linear space and 

t'' : X ^ L^{n) 

linear operators. A sequence {u^} in L^(f2) is said to two-scale converge to 
Uq ^ L^{n X Y) with respect to {t''} if, for any v G X , 

\im J u'^{x){t'^v){x) dx ~ J J uo{x,y)v{x,y) dy dx. 
n n Y 

In order to achieve a compactness result like Theorem[31 we need Dcfinition[5l 

Definition 5. Assume that Y is an open bounded subset of M.^^ . Let X C 
L^(il X Y) be a normed space and 



t'' : X ^ L^{n) 
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linear operators. Then {r''} is two-scale compatible with respect to X if there 
exists C > independent of h such that, for any v Cz X , 

lim ||T''i;||i2(o) s$ C||i;||L2(nxy) 

h — >oo 

and 

I|tS||l2(o) ^ 

Furthermore, X is called admissible with respect to {t'' }. 

Using Definition [5l wc have the compactness result according to Theo- 
rem [6] [10. 

Theorem 6. Assume that Y is an open bounded subset o/R^^. Let {r''} be two- 
scale compatible with respect to X , a separable Banach space dense in L^(ilx Y). 
Then every bounded sequence {u'*} in L'^{il) has a subsequence that two-scale 
converges with respect to {t'*}. 

Using weak convergences, we can introduce a stronger version of two-scale 
compatibility than presented in Definition [51 see Definition [7] 

Definition 7. Assume that Y is an open bounded subset o/R*^. Let {t'*} be 
two-scale compatible with respect to the linear space X C L^(fl x Y), and let 
{u^} be a bounded sequence in L'^{n) two-scale converging to Uq. Then {t^} is 
strongly two-scale compatible if, for any v G X , 



Y 

and 



UQ{-,y)dy inL'^{n) 



T^v v{-,y)dy inL^in). 



Compare Definition [3 with the result of Proposition [2 

An alternative generalization of two-scale convergence is scale convergence 
introduced by Maria Luisa Mascarenhas and Anca-Maria Toader in 2001 [5]. 
The definition of scale convergence is according to Definition [8] 

Definition 8. Assume that Y is a metrizable compact space, /i a Young measure 
onilxY , and L^^iVl x Y) is the space of all functions with fi-integrable square. 
Furthermore, let {a^} be a sequence of fi-measurable functions a'^ : ^ Y . A 
sequence {u''} in L^{n) is said to scale converge to uq E x Y) with respect 

to {a''} if for any v G L'^{n; C{Y)), 

lim / u^{x)v(x,a'^{x)) dx = / UQ{x,y)v{x,y) d^(x,y). 

h^oo J \ ' / J 

n nxY 
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A special case of scale convergence is achieved by choosing the same mea- 
sure and same class of admissible test functions as in two-scale convergence, 
namely the Lebesgue measure A and Cper(^)), respectively. This leads 

to Definition [5] [131 Definition 27] (see also [51 Definition 11]). 

Definition 9. Assume that Y is the unit cube in M*^. Furthermore, let {a^} 
be a sequence of functions a'^ : fl ~* Y . A sequence {u''} in L^(f2) is said 
to X-scale converge to uq G L'^{fl x Y) with respect to {a'*} if, for any v € 

^lim j {x)v {x , a'^ [x)) Ax — j j uo{x,y)v{x,y) dy dx. 
n n Y 

A thorough treatment of how {a'*} could be chosen to obtain strong two-scale 
compatibility is found in the doctoral thesis [T31 Subsection 2.4.2] of Jeanette 
Silfver (see also F5]). We reproduce the results below. 

Following [13i Subsection 2.4.2], we omit cases where M ^ N letting 

a'' : ^ 

be a continuous bijection. Furthermore, {Y^}°^i is a covering of with unit 
cubes, and {Yj!}^^^ a covering of Y^ with cubes with side lengths ^ in such 
a way that Yj! are F-periodic repetitions of a cube Yk C Y. Assuming that 
r2 C has a Lipschitz boundary, we define 

n'; = {a'')-\Y^)nn. 



We also assume, for each h, that there exists a finite set q{h) C Z+ such that 



= U,g„(/i)r2''. Similarly, introduce 



Finally, it is assumed that flj C Nj.{h){xhj) for some ball Nr{h){xh,j), centered 
at some Xhj G Vl^, with radius r{h) — > as /i — > cxd. Given the setting above. 
Definition [TUl makes sense [131 Definition 28] (see also [SI Definition 12]). 

Definition 10. Suppose that for all cubes Yk 'ZY and any i^j, j £ q{h), such 
that {a'^)^^ (Y^) n fl does not intersect dft, 



<e (1) 



where e = e{h) — > for h ^ oo. Then {a^} is said to be asymptotically 
uniformly distributed on Q. 

This is sufficient to obtain the strong two-scale compatibility as promised, 
see Proposition [TT] [T3l Proposition 30] (sec also [5, Proposition 15]). 
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Proposition 11. Assume that f2 has a Lipschitz continuous boundary. Let 
{a^} be asymptotically uniformly distributed on fl. Then {r''} defined by 

is strongly two-scale compatible with respect to L^(jl; Cpcr{Y)^. 

Note here that the admissible space in L'^{QxY) is in this case L^{^; Cpcr(y)) ■ 
As a consequence of PropositionfTTl we have Corollary [T^ 13, Corollary 31] (see 
also Corollary 16]). 

Corollary 12. Assume that 51 has a Lipschitz continuous boundary. Let {a^} 
be asymptotically uniformly distributed on il, and let {w''} strongly converge to 
u in L'^iyi). Then, up to a subsequence, X-scale converges to u with respect 
to {a^}. 

Note that the second scale dependence of the A-scale limit vanishes, just like 
how it works for two-scale limits in the case of strong convergence [7] . 

3 Convergence for sequences of A/^-tuples of par- 
tial derivatives 

Since sequences of partial differential equations in the context of homogenization 
typically involve A^-tuples of partial derivatives of solutions, i.e., Vu'', we need 
to investigate the two-scale limits for these. Indeed, for traditional periodic 
two-scale convergence, we have Proposition [T3] [Tl 110] . 

Proposition 13. Let {«''} be a bounded sequence in H^{Q,) such that the 
strong limit in L'^{fl) is u. Then, up to a subsequence, there exists Ui G 
L^{n; Wp„{Y)) such that for any v G L^{n; Cp^,{Y)^), 

N 



f du 
lim / y — — {x) Vi{x,hx) dx 

h^ooj OXi 



n Y 



4=1 



It should be noted here that bounded sequences in H^{^) strongly converge 
in L^(ri), and we know that strongly convergent sequences in L^(ll) have a 
subsequence with a two-scale limit with vanishing second scale. 

A deciding step towards the homogenization of certain non-periodic prob- 
lems is to prove the corresponding result for A-scale convergence. We have 
Proposition [14] proved by Jeanette Silfver in [13l Proposition 35] (see also [3 
Section 4]). 
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Proposition 14. Let X be a Banach space for which 
and {a^} a sequence of functions 



which are continuous and bijective. We assume also that {t''} defined by 



r V = v\ • ,a 



is strongly two-scale compatible with respect to X . Furthermore, let Z C L^{^', Lp^^{Y)^ 
be a Banach space such that Z D Vift; C^^^{Y)^) is dense in Z , and let Z 
be its orthogonal complement in L^iyi] L'^^^^{Y)^^ . Assume also that, for any 
vGZnV{n;C^,,{Y)^), 

N N ^ cj h 

EES:^(-"''(-))^-0 ^nL\n). (2) 



Then, for any bounded sequence {u'^} in H^{il) there exists a subsequence 
such that ^ u in L'^{il) and, for any v ^ X , 

\im J u^{x)v(^x,a'''{x)) dx = J J u{x)v{x,y) dy dx, (3) 
n n Y 

and there exists Wi G Z-^ such that, for any v £ X^ , 

N 



lim / — — (x)vi(x,a'^(x)) dx 

h^ao J ^ dXi ^ ' 



N 



E ( 'd^^^^ +wi,tix,y) jvt{x,y)dydx. (4) 



n Y 



Proof. Property ([3]) is given by Corollary [121 By using Green's formula twice, 
using assumption ([2|), and utilizing density, property ([4]) follows from orthogo- 
nality. □ 

Remark 15. The original proof of Jeanette Silfver with all details is found 
in IM Subsection 24.2]. 



4 Elliptic homogenization 

In this section we present a homogenization result for elliptic problems governed 
by the sequence {a'*} introduced earlier in this paper. Such results where first 
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published in [T3] , and in Sections \E\ and [H we present a special case in two 
dimensions and perform a numerical experiment, respectively. 

The foundation upon which modern homogenization rests was erected in 
1967-68 when Sergio Spagnolo developed the concept of G-convergence of se- 
quences of symmetric matrices [131 [TS], which has later been generalized by 
Frangois Murat to H-convergence where the symmetry assumption is dropped, 
but at the cost of an imposed requirement on the sequence of flows [9] . We begin 
by introducing a space of matrix valued functions according to Definition 1161 

Definition 16. Ifoo>r^s>0 and O C is open, the matrix valued func- 
tion M e ig said to belong to A4{r,s; O) if the system of structural 
conditions 



sup 



N 



(hounded) 



N N 



inf (.ifTiij^j 5^ (coercive) 



1=1 3 = 1 



is satisfied a.e. in O. Furthermore, the space Mpci{r,s; Y) consists of those 
functions in A4(r,s; which are Y -periodic. 

We can now give Definition [T71 defining the notion of H-convergence. 

Definition 17. Let {A^} be a sequence in A4{r, s; il), and let B ^ A4{r' , s' ; fl). 
Furthermore, assume that, for any f £ H^^{n), the sequence of solutions {u'*} 
in HQ{fl) to the sequence of problems 



N N 



^ ^ ^X^ 

1=1 3 = 1 



in il 
on dn 



(5) 



satisfies 



N 

E 



dx. 



u in Hq{Q) 



N 



du 

3 = 1 J 



in L'^ifl) 



N 



(6) 



where u is the unique solution to 

N N 



EE 



1=1 3 = 1 



dxi 



dx, ' " ^ 



in U, 
on dfl 



(7) 



Then {A'^} H-converges to B. 
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Remark 18. In the literature, H-convergence is often called G- convergence. 
In this paper, we clearly separate the general notion of H-convergence from the 
special case of G-convergence which deals with symmetric matrices exclusively. 

When an H-limit has been found, one has homogenized the sequence of 
partial differential equations (O, and ([7]) is the homogenized problem. For 
sequences of matrices which are periodic, we get Theorem [19] for the homoge- 
nization in this case [TOl [1] . 



Theorem 19. Suppose A E A4pcr{r,s; Y). Then {A{h-)} H-converges to B 
given by 



k=l ^ " ' / i,j = l,...,N 



where z € Wpcr(5^)^ uniquely solves the local problem 

N N 

rl I / \ i \ 

= inY. (8) 

j=l,...,N 

We note that it is the term containing the local problem solution z which 
makes the H-limit B to deviate from the average of A over each periodicity cell 
Y. It is due to the ui term in the two-scale convergence for partial derivatives 
and is obtained through a simple separation of variables of ui with Vu as the 
a;-dependent part and z as the y-dependent part. 

In the remainder of this paper, A^ in the sequence of problems (O will be 
the composition Ao . 

Below we present the result on non-periodic homogenization from [13j an- 
nounced in the beginning of this section. We start by giving Definition [20] [13|, 
Definition 51]. 

Definition 20. The sequence {a'*} is of type if the following three condi- 
tions hold: 

(i) The sequence {t^} of operators : X ^ L'^{fl) defined by 

T''v = v{-,a''{-)), veX 
is strongly two-scale compatible with respect to a Banach space X for which 
V{n; Cp- (r)) C X C L%,iY)). 

(ii) There exists a sequence {p^} strongly convergent to zero in H^iyi) such 
that 

"'(S).,,.. »-('">-.^- ~ '»^'("'"^" 

where 11 is a diagonal matrix with 



li=l,...,N 



(CO.=i,... V e L<-{nf. 
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(in) There exists a Banach space Z C for which Z n 'Difl; C^^(y)^) 
dense in Z , and such that, for any v £ Z C] C^^.{Y)^'j , 



From Definition [20] and Proposition I14| we have the homogenization result 
of Theorem ED [131 Theorem 52]. 

Theorem 21. Assume that has a Lipschitz continuous boundary and that 

AeCp„.(r)^><^nA^(r,s; R^). 
Let {a^} be of type H^2(q.(;; {Y})' ^^^^ {Aoa'^} H-converges to B given by 

N ^ X / „ N 



+ ,y) j dy \ 



u G Hq{^) being the weak limit of the sequence {u^} of solutions to ([5]), if 
u e H(){^) and wi € Z-^ C L/^^^iY)'^^ uniquely solve the homogenized 

problem 



\ / EE )^(^)^2/dx = / /(.T)u(x)da; (9) 



for all V e iJQ(f2), and, for each x eO,, the local problem 

N „ N 



for all V e WpcriY)- 

Proof. First we make a weak formulation of the sequence of problems ([5]) . 

The homogenized problem ^ is obtained by passing the weak formulation 
to the limit, using Proposition [HI 

The local problem (|10p is derived by first choosing appropriate test functions, 



then utihzing the fact that {a'*} is of type H^2(q. ^-^^^.^y)), and finally employing 
the variational lemma. □ 

Remark 22. The original proof of Jeanette Silfver with all details is found 
in pJl Subsection 3.5.1]. 
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5 A non-periodic and two-dimensional example 



In order to justify the concept of homogenization of A-structures, we take a look 
at a non-periodic two-dimensional example (thus, fixing N = 2 from now on). 
Let {a'^} and be given by 

a^ix) — hxi „ 
Qfj (a;) — hx2\x2\ 

and f] — (ai, 6i) x (02, 62), respectively, where 6.^ > a.^ > 0, i = 1, 2. Thus, H is 
an open interval in the first quadrant of M?. In Figure [T] we depict in the first 




Figure 1: Some entry of A o a'', /i = 1, in the first quadrant. 

quadrant the behaviour of some entry of Ao for some (0, l)^-periodic matrix 
A. (To be specific, we have chosen an entry on the form 1-1- 1 siuTryi sin7r2/2|^^^-) 
Note how the oscillation frequency increases with growing X2 , while the period- 
icity is preserved in the xi direction, as expected. 

To homogenize the sequence {A o a'^}, we must first check that {a''} is 
asymptotically uniformly distributed on fl. Indeed, we have Proposition 1231 

Proposition 23. The sequence {a''} given in (jlip is asymptotically uniformly 
distributed on 51. 

Proof. We define = j+Y and F^, = i(fc + F), where j e Z^, n G Z+ and 
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fc G {0, . . . ,n — 1}^. Furthermore, = j + Yj^. Clearly, there exist bijections 
(j: {1,2,...} 

\/C: {l,...,7i2}^{0,...,n-l}2 

such that {Y-^''^^}°°^^ is a covering of R'^, and {^^If)'' })!_]^ is a covering of 
Y'^^^K This merely shows that the way of enumerating the squares in this proof 
is equivalent to the way in Definition (TUl 

We have, for each j — (ji, j2) G 7? and fc = (fci, G {0, . . . , n — 1}'^, 



and 



where 




ft'* 



and 




fci+i 




In order to cover O, we can for each fixed h apparently do this with a finite 
union of fij sets. The smallest ji required for the covering will go like 0{h) and 
the smallest j2 like 0[\/h). For large h, the diameter of the fij sets will go 
uniformly like C'(--^), so we can fit them in balls of equal radius which tends 

to zero as ft, ^ oo. We realize that what is left to check is ([T]). 

For covering sets, let j be chosen such that {a'^)~^{Y^) does not intersect 
dV,, i.e., 

For such j's we have the Lebesgue measures 



[Vh + 1 - VT2) 
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and 



1 1 

hVhn 



J2 



n 



Thus, 



J2 + 



k2 + l 



J2 + — 

n 



y/j2 + 1 - V32 
Vh + 1 + Vj2 



1 



]2 



J2 



Since the smallest j2 goes like O(V)i), we must uniformly have that 



for any k £ {0, 



, — 1}^, a,s h ^ oo. This implies 



< e, 



where e = e(/i) ^ as /i — > oo. 

In virtue of Proposition [TTl we thus get Proposition [24l 
Proposition 24. The sequence {t'*} defined by 



T V — v\ • ,a 



□ 



where {a^} is given by (jlip . is strongly two-scale compatible with respect to 

L\n- Cper(y)). 

Proof. Use Proposition [23l together with Proposition [11] □ 
Define the 2-tuple C according to 

fCi(a:) = l 

\C2W =2X2 ' 

and define the diagonal matrix 11 element-wise by 

(^"W).=l,2 = (C^(^))^=l,2- 

Proposition [25l below shows that {a^} is of right type to use in Theorem [2l1 for 
the homogenization to work out properly. 



(12) 
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Proposition 25. The sequence {a''} given by is of type H^2(q.(7 (Y))' 

where ^ is given by (I12p . 

Proof. We must check conditions (i)"(iii) of Definition [20l Condition (i) is an 
immediate consequence of Proposition 1241 For condition (ii), we note tliat 

- — [x) = 2/1x2 , 
0x2 



, dx 



and by choosing 



for which p'* ^ in iJ^(ri), we get 



We obviously have that C G L°°(ri)^. FinaUy, we must check condition (iii) 
of Definition 1201 In this context, it means that we must find a Banach space 
Z C L'^{n- Cpc,{Yf) for which Z n V{fl- C^^^{Yf) is dense in Z, and such 
that, for any w G Z n V{VL- C^^^{Yf), 



dvi 



'.i:i(-""(-))c, 

as /i — > 00. If we define 



in L^{n), 



(13) 



(14) 



we clearly have, for any u e Z n 'D(il; C^^{Y)'^) , a satisfied weak convergence 
(fT3l) . What is left to prove is the density oi ZnV^D,; C^^^{Y)'^) in Z, where Z 
is given by (I14p . We apparently have 



(15) 



and since V(n; C^^^.{Yf) is dense in L'^(n; L^^^iYf), it is easy to check that 
we have the desired density of Z n X>(fi; C^^AY)'^) in Z. □ 

Remark 26. H^e /lave silently used the fact that the vanishing divergence in 
the definition for Z and the expression for Z D 'Dijl; 6*^^(1")^) in (jl4[) and 
(jl5[) . respectively, does not upset the inheritance of the density property for 



v[n-c^,,{YY) inL'[n- Li,,{Yy 
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To obtain the orthogonal complement, we need Lemma [771 [TU] , 

Lemma 27. Assume that Y is a unit cube in M.^ . Let f e Lp^j.{Y)^ be 
orthogonal to the space 



{9 e C-AVr : E S - 0} 



of divergence free functions. Then, for some h G Wpcr(l"), 



^ V%/i=l,...,7v' 

We are now ready to characterize the orthogonal complement. 

Proposition 28. The orthogonal complement C L/^^^iy^^^ of Z, 

defined in the proof of Proposition \25\ is 

Proof. Suppose v G Z and wi G Z^. Then, by definition, 

2 

0==y" J ^ v^{x,y)wi^i{x,y)dydx 
n Y '=1 

2 



n Y 

Since 



Lemma [27l implies that, for some ui{x, • ) S VFpcr(i^) a.e. x E ft, 
Hence, for some ui{x, • ) € Wpcr(i^) a.e. a; € 17, 

We know that Z^ C ^^(j]; ^^^.(y)^), so wi G ^^^.(F)^) holds. This 

implies 

< max{l, 2^}|lwi|li2(f^.i2_^^(y)2) 

< CX), 
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where we have recalled that = (ai, 61) x (02, 62)- Thus, ui G L'^{^; Wpcr{Y)), 
and we conclude that 



and we are done. □ 

We can now formulate a preliminary homogenization result in Proposi- 
tion [21 

Proposition 29. Assume 

AeCpcr{Yf''^nM{r,s; R^) 

and let {a^} and C be given by (llip and (jl2p . respectively. Then {A o a^} 
H- converges to B given by 

J = l ^ \y i = l / i=l,2 

u G i?o(f2) being the weak limit of the sequence {u^} of solutions to (O, if 
u G H}^(yi) and ui G IVpcr(i^)) uniquely solve the homogenized problem 



n Y 



f{x)v(x) dx 



for all V G Hq{Q?), and, for each x G 51, the local problem 



/^^a,,(y)(^(.)+O(a:)|^(x,y))c.(x)^(y)dy = (17) 
/or aZ/w G VKpor(i^). 

Proof. This is an immediate consequence of Theorem [21] together with Propo- 
sition [211 and Proposition [551 CH 

It is possible to improve Proposition [2H| to yield an explicit homogenized 
matrix and a local problem of the same type as in Theorem 1 191 Indeed, we have 
Theorem [301 

Theorem 30. Assume 

AG Cpcr(i^)^''^n Al(r,s; R^) 
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and let {a^} and ( be given by ([TT|) and (fT^ . respectively. Then {A o a''} 
H- converges to B given by 

iK), I jl^^k{y)(5u,+(:k^{-,y))d^ , (18) 



where z E L°° (jl; T4^pcr(^)^) uniquely solves, for each x G ft, the local problem 
- (EEC'(2:)^(a.fe(4,+aw|g(a;, O))^ =0 inY. (19) 



Proof. The proof will be performed in four steps, where the first step introduces 
an ansatz, the second and third steps derive the homogenized matrix (|18[) and 
the local problem (fT9|) . respectively. In the last step we prove the uniqueness of 
the solution to the local problem. 

Step (i): Ansatz. Let u G -ffo(f^) and ui G L^{fl; W^por(i^)) solve the system 
of equations in Proposition I29i and make the ansatz 

du . , 



^3 



where z G L°° {p,\ Wpcr(y)^)- We will see later why it is necessary that we must 
constrain ourselves to z( • , y) e L°°{flY a.e. y G M^. 

Step (ii): Homogenized matrix. Let us first derive the expression (jlSp for 
the homogenized matrix. From and ([^0]) we get 

^J = l ' = l ^ fe=l ^ / i=i^2 

which is satisfied if 

(^'^).,,=i,2 ^ (/E «»fc(y)('5fe. +a|g(-,y)) 

where we note that B e L°°(ri)^^^, which requires z{- ,y) G L°°(ri)^ a.e. y G 

S'tep (^Mij.' Local problem. Next, let us derive the local problem (|19p . Fix 
some V G VFpcr(5^) to be used in the local problem ([TT]) in Proposition [211 whose 
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left-hand side becomes, for each x G fl, 



E E (y) [ + (^) ^ E (^)^^- y) J (^) (y) 

7 — 1 -i^ i— 1 fc— 1 



This must be zero, which is the case if, for each x e 17, 

[jilil^^ (^) ^ + a (x) || (x, y)) dy J =0. 

\y i—^ k—1 / j—\^2 

By partial integrating and using the divergence theorem, we obtain, for each 



kj + Q^{x)^{x,v))dS\ 

\QY ^^=1 ' ' j = l,2 

- (/ «(y)EEC.(^)^(a.fc(y)(4,+a(:^)|^(2:,y))]dyJ =0, 

\y i—\ k—1 / j — 

where n is the unit outward normal to dY. Since v is y-periodic, the surface 
integral vanishes, and we are left with, for each x G ft, 

-(/ %)EEC.(x)^(«./c(y)(4,+a(^)|^(^,y)))dy) =0, 

which certainly is satisfied if, for each x € Q, 

-(i:i:c.(-)|-(-(*«+&w|£)])^^_^=». 

and we have shown (fT9l). 



J = l,2 



Step (iv): Uniqueness. It remains to prove the uniqueness of the solution to 
the local problem. Fixing x £ fl, we can define a new, rescaled, y-variable y?(^) 
by letting 

The new local problem can be written 
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effectively a "classical" local problem, where the 2x2 matrix A'^^^^ and the 
2-tuple z^(^) are given according to 

\ z^t^) o y';(^) = z{x, .) ' 
and y?(^) = (0, 1) X (0, Since 

uniqueness is ensured just as in the "classical" case. □ 

Remark 31. In Step (iv), the fact that Y^'^^\ for each x £ Q, is a rectangle 
rather than the unit cube will, of course, not spoil our argumentation. 




0.5 1 1.5 2 



Figure 2: The scalar factor of ^ o a'', h = 3. 

6 A numerical illustration 

As an illustration of the theoretical results of Section [51 consider the sequence of 
problems ([5]) with = Ao where A is given as a product between a scalar 
factor and a unit matrix according to 

i'^viy)) ^J=l^2 = + ^sin27rj/i sin 27r2/2) (^»j)»j=i,2, y e M^ 

and let = {S, 2)^ where S >Q. Apparently, 
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Ql ^ ^ ^ 1 

0.5 1 1.5 2 

Figure 3: The diagonal entries of B. 

so Theorem 1301 is applicable. Furthermore, 

((fly o a'')(x))^^^.^^ 2 = (l + ^sin27r/ia;i sin27r/ia;^)(%)zj=i,2, xG ((5,2)^ 
see Figure [21 By using (fTg]) . the local problem for z becomes, for each x £ fl, 

- (EE = inr, (21) 

where the matrix C and the 2-tuple g are given by 

cii{x,y) = (l + ^ sin27r?/i sin27r?/2) 
C22(a;, y) = (l + ^ sin27ryi sin27r?/2) 

and 

gi{x,y) = fvr cos27r?/i sin27ry2 
g2{x,y) = ^7ra;2 sin27ryi cos27ri/2 ' 

respectively. Solving (j2ip numerically (effectively we have a one-parameter fam- 
ily, with respect to X2 , of partial differential equations to solve) and then com- 
puting the homogenized matrix through HH]), we get that B is diagonal with 
non- vanishing entries, functions with respect to X2 only, given according to 
Figure [3l Note the interesting feature that 6ii|a;9=i/2 — ^221x2=1/21 where B 
obviously is proportional to the unit matrix, i.e., along the line X2 — ^, the 
homogenized matrix is isotropic. The heuristic explanation to this is simple; 
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for large h the mapped periodicity cells in the vicinity of the line X2 = 2 s-re 
near-perfect squares. Of course, it is crucial that A is isotropic to begin with in 
order for the map Ao to exhibit a near-isotropy property on such mapped, 
near-perfect squares. 
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